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ABSTRACT 

We find an exact solution of Kerr-Newman-de Sitter type on the braneworld(4D) of the 
DGP model. When a constant 4D Ricci scalar is assumed, only zero(flat) and a positive(de- 
Sitter) values satisfy the Hamiltonian constraint equation coming from the extra dimension. 
With a Z 2 -symmetry across the brane and a stationary and axisymmetric metric ansatz on 
the brane, we solve the constraint equation exactly in the Kerr-Schild form with de-Sitter 
background. In the de-Sitter background this Kerr-Schild solution is well behaved under 
Boyer-Lindquist transformation: the constraint equation is preserved under the transforma- 
tion and so is the solution. In the non-rotating limit we show that this Kerr-Newman-de 
Sitter solution has the characteristic of accelerated expansion of the braneworld universe. 
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1 Introduction 



Recent astronomical observations of Type la supernovae suggested that our universe is ex- 
panding at an accelerating pace pQ. In parallel with this development, the idea that our 
universe may be a brane embedded in some higher dimensional space also became quite 
popular [2]. One of the models along this line proposed by Dvali, Gabadadze and Porrati 
(DGP) [3] is known to contain a branch of solutions exhibiting self-accelerated expansion of 
the universe [I] . In this model, the acceleration takes place without having the cosmological 
constant on the brane. Although there have been many works and much interest in the DGP 
model, its exact solutions are not much known so far: It's approximate Schwarzschild solu- 
tions had been obtained in [H El El [HI E] , and an exact Schwarzschild solution on the brane 
was obtained in pi)]. Recently, an exact solution for charged black holes on the brane was 
obtained in [11], and we worked out a particular solution for charged rotating black holes on 
the brane in [12]. However, in our previous work the metric in our modified Boyer-Lindquist 
coordinates could not be fully diagonalized: it contains t-9 and Q-ep cross terms besides the 
usual t-cf) cross term for rotating solutions. This causes a difficulty in interpretation of the 
structure of the horizon. 

On the other hand in [12] we also noticed that the particular solution we obtained shows 
the characteristic of Kerr-Newman-de Sitter black hole in general relativity. Based on this 
observation, we further investigate along this direction and in this paper we obtain an exact 
solution of Kerr-Newman-de-Sitter type with the metric in conventional diagonalization with 
t-(p cross term only. 

In doing this, we first notice that only two solutions, flat and de-Sitter type geometries, 
are allowed by the constraint equation if constant curvature scalar is assumed. Then we 
solve the constraint equation in the Kerr-Schild form [13] with de-Sitter background and 
find an exact solution. This works well with the Boyer-Lindquist type coordinates [H] in 
the conventional sense in which only t-cf) cross term appears in the metric. This solution 
also matches well for a cosmological solution of accelerated expansion: Following the work 
of [15], we show that our de-Sitter type solution exhibits the characteristic of accelerated 
expansion in the non-rotating limit. 

This paper is organized as follows. In section 2, we set the action and equations of motion 
of the DGP model following the approach of [16J. In section 3, we obtain a Kerr-Newman 
type exact solution on the brane with de-Sitter background in the Kerr-Schild form and 
transform it to the Boyer-Lindquist coordinates. In section 4, we show that our de-Sitter 
type solution matches a cosmological solution with accelerated expansion in the non-rotating 
limit. We conclude in section 5. 
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2 Field equations on the brane 



The DGP gravitational action in the presence of sources takes the form [3] 

S = Ml J d 5 xV=9 {5) R + J d A xV^h (M 2 P R + L matter ) , (1) 

where R and ^R are the 4D and 5D Ricci scalars, respectively and L matter is the Lagrangian 
of the matter fields trapped on the brane. Here, the (4 + 1) coordinates are x A = (x^, y{= 
x 5 )), fi — 0,1, 2, 3, and g is the determinant of the five- dimensional metric g AB , while h is 
the determinant of the four- dimensional metric h^ v = g^ix^^y = 0). A cross-over scale is 
defined by r c = m" 1 = Mp/2M^. There is a boundary(a brane) at y = and Z 2 symmetry 
across the boundary is assumed. The field equations derived from the action (TjQ) have the 
form 

^G AB = ^R AB - \g AB ®R = 4\ - ( X ab + T AB ) 5(y), (2) 



^g 

where n\ = Mp 2 and = M~ 3 , while X AB = —5^5 b G^ u /kI and T AB = 5^5 B T^ U is the 
energy-momentum tensor in the braneworld. 

Now, we consider the metric in the following form [Toj IT7] , 

ds 2 = g AB dx A dx B = g^{x, y)dx^dx u + 2N fl dx>*dy + (N 2 + g^N^N^dy 2 . (3) 

The (/i5), (55) components of the field equations (j2J) are called as the momentum and 
Hamiltonian constraint equations, respectively, and are given by [Tol IT?] 

V V K^ - V,K = 0, (4) 

R-K 2 + K^K^ = 0, (5) 
where is the extrinsic curvature tensor defined by 

K»» = ^{d y g,u - V M iV, - V V N„), (6) 

and V M is the covariant derivative operator associated with the metric g^. 

Integrating both sides of the field equation ([2]) along the y direction and taking the limit 
of y = on the both sides of the brane we get the Israel's junction condition [18] on the Z 2 
symmetric brane as follows [TTj : 

G>„ = k\T^ v + m c (K^ u - h^K). (7) 
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In this paper, we take the electro-magnetic field as the matter source on the brane. 
Plug the junction condition ((7|) back into the constraint equations (jlj) and we find that 
the momentum constraint equation is satisfied identically, while the Hamiltonian constraint 
equation is given by 

R^R^ - ^ 2 + m 2 c R + <T^ V - 2k\R^ v = 0, (8) 

where we used T = = 0. Here we would like to note that when the spacetime has constant 
curvature scalar only two values are allowed, R = 0(flat) and R = 12m 2 c (de-Sitter). This 
can be seen as follows: For R = 4a with some constant a, one can set R^ := ah^ v + k^T^. 
Plug this relation back to the constraint equation (jHJ), one gets a = or a = 3m 2 .. 
Finally, combining the Einstein equations in the bulk(|/ ^ 0), 

^G AB =^ R AB - X -g AB ®R = 0, (9) 
with ((7|) we arrive at the gravitational field equations on the brane [11] 



k\ 1 
Gfiv E^ v —(T^Tpi/ "^h^i/TpuT^ ) 



1 2 1 2 1 
2(R P fiRpv — -^RRfiu + -rh^R — -hp jU R pa R pa ) 

+^L(RP^T pi/ + T^R pu - hlT^ - h^R^T^), (10) 

where is the traceless "electric part" of the 5-dimensional Weyl tensor ^Cabcd [19J and 
m^ 1 = k\/2k\. In what follows we shall set k\ = 8it. 

Here, we would like to check the existence of the gravitational effect from the extra 
dimension. Use the fact that the Hamiltonian constraint (jHJ) is satisfied by the relation 

Rixu = oih,„, + itfFf,,,, (11) 

for a = or 3m 2 , then the gravitational field equations fllOp on the brane reduce to the 
simple form 

R pu = ah pu — E pu . (12) 



Comparing (1121) with (TTTI) we can identify the projected Weyl tensor E pv with the energy- 

E^u = — k 4 T^v (13) 



momentum tensor T pu on the brane as 
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This shows that for the two special cases of R = and R = Ylm^ the tensor E pv has only 
contribution from the energy-momentum tensor on the brane. Therefore our field equations 
on the brane, (1101) contains no effect due to the extra dimension when the spacetime has 
constant curvature scalar. This is different from the RS model case [20] , where the so-called 
tidal charge effect due to the extra dimension appears. 



3 Kerr-Newman solutions on the brane 

In this section, we consider the case in which the brane contains a Maxwell field with an 
electric charge. In general, rotating black holes on the brane can carry electric charges. We 
assume that the Maxwell field on the brane is described by a solution of source-free Maxwell 
equations. In that case, the trace of the energy-momentum tensor for the Maxwell field on 
the brane vanishes, and we have to solve equation ^ and the Maxwell equations: 

g^D,F ua = 0, (14) 

D^F V(j] = 0, (15) 

where is the covariant derivative operator associated with the brane metric h^ v . However, 
since ( f!5|) is satisfied identically once the field strength is expressed in terms of the potential 
one-form, we only have to solve equations ([8]) and (JHJ). 

To solve equations ([H]) and (|T4^ . we start with a stationary and axisymmetric metric 
ansatz describing a charged rotating black hole localized on the 3-brane. We write the 
metric in the Kerr-Schild form [13] with the metric expressed in a linear approximation 
around the background: 

9w = 9pu + vk^ky (16) 

where v is an arbitrary scalar function and k^ is null, geodesic with respect to both the 
background metric g^ u and the full metric g^ u . 

In general, it can be shown that if the vector field k^ is null-geodesic, then Ricci tensor 
of g^ u is related to that of g^ v by [21] 

R" u = R\ - s%W v + l -D p D vS w + \lPTP8 vp - \D p D p s^ (17) 

where = vk^k u and is the covariant derivative with respect to the background metric 
g^. Thus, in our case if g pv satisfies R pv = Ag^, then the full metric g^ u satisfies the 
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Einstein equations R^ u = Ag^ u + k\T^ v with the same cosmo logical constant, provided that 
s^ u satisfies the linearized Einstein equations with respect to the background metric g^ u . 

In this paper we will mainly deal with spacetime(4D) with constant scalar curvature. 
In the previous section we saw that this choice leaves us two possibilities, flat or de-Sitter 
spacetime. In the flat (R = 0) case, the solution for equations (jHJ) and (JHJ) is given by the 
usual Kerr-Newman solution in general relativity [12]. Thus from now on we concentrate 
ourselves on the de-Sitter case, and choose the de-Sitter background metric in the Kerr-Schild 
form to solve equations (jHJ) and (fill) . 

Introducing the coordinates x M = {r, r, 9, ip}, we write the metric in the Kerr-Schild form 
as li- 



ds 2 



where 



-(1 - Ar 
+H(r,9) 



-dr 2 + 



^dr 



(1 -Xr 2 )(r 2 + a 2 ) 



, 2 £ Jn2 (r 2 + a 2 ) sin 2 6 , 2 
dr 2 + —d9 2 + - ^ dip 2 



(1 - Ar 2 )(r 2 + a 2 ) 



a sin 2 
dr — — dip 



;is) 



E 



2 , 2 2 n 

r + a cos y, 
1 + Xa 2 cos 2 0, 
1 + Xa 2 . 



(19) 



Here a is related to the angular momentum of the black hole and A is a constant parameter 
determining the curvature of the background metric. 

With the de-Sitter metric (TTBl . the Maxwell equation (fill) is satisfied if the potential 
one-form is given by 



Qr 

A^dx^ = —^(—AgdT + asm 2 9dip), 



(20) 



where the parameter Q is the electric charge of the black hole. The nonvanishing components 
of the electromagnetic field tensor F^ v for the above potential one-form A is calculated to be 



_ QA e (r 2 -a 2 cos 2 fl) _(1_- Xr 2 )Qra sin 29 



5£ 2 

(r 2 - a 2 cos 2 9)Qa sin 2 9 



5£ 2 

-2 i „2 



r 2 + a 2 )Qrasin29 



rip 



^YJ2 ' •* v<p ^Y? 

Using (El]) and setting the scalar function H in the metric ansatz (1181) as 



(21) 



H(r, 



rP{r,9)/E, 



(22) 
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the Hamiltonian constraint © can be expressed as 







72A(A - m 2 )E 2 + 6(m 2 - 2A)E(2Z + rZ r ) + r 2 Z 2 

4 



S/' I 2 fl2 COs2 $ 



E 
12Q 2 



ZZ r . + — (r 4 - 4a 2 r 2 cos 2 + a 4 cos 4 fl)Z 2 
E 2 



[2(r 2 -a 2 cos 2 #)Z-rEZ r ] + 



24Q 4 



(23) 



E 2 L_v " y ~ '~"' J ' E 2 ' 

where Z = dP/dr and Z r = d 2 P/dr 2 . For A = and a = 0, the above equation reduces to 
the equation for charged black holes in a flat background [TT] . 
In order to find a general solution we solve ( 1231) for Z r : 

6Q 2 /4 + 2a-2a 2 * 



— rZ r 



r 2 (l +a) 



1 + 2a + a 2 



Z + 3(m 2 c -2A)r 2 (l + a) 



±3m 2 r 2 (l + a)\/l + Au + ^w 2 , 



(24) 



where a(r, 6) = a 2 cos 2 6>/r 2 and u = [Q 2 — r 2 (Z + 3Ar 2 (l + a) 2 )]/m 2 r 4 (l + a) 2 . Introducing 
a new function U by putting Z as 



Q 2 /r 2 - 3m 2 r 2 (l + a) 2 U(z, 6)/2 - 3Ar 2 (l + a)' 



(25) 



where r = roe 2 , we express ( 1231) in terms of {7 
"2Q 2 3m 2 r 2 



+ 



(2(1 - a 2 )U + (1 + a) 2 t/ 2 ) + 6Ar 2 (l - a' 



20 2 3 

+ -m 2 r 2 (4 + 2a - 2a 2 )U + 6Ar 2 (l - a 2 ) + 3m 2 r 2 (l + a)(l ± /) 



where / = a/1 + 6C/ + 3£7 2 . Solving f[2"6"]) in terms of U z = dU/dz, we get 

2 



17, 



1 + a 



[l + U±f). 



, (26) 



(27) 



Here we would like to mention two special solutions, U — and U = —2 for minus and 
plus signs, respectively, for the last term in the above equation. Integrating ( 1271) using the 
relation r = roe 2 , the solutions for U can be given implicitly by 

k S 4 = (1 + 317 + /) 

1 £/ 2 (3 + 3£/ + A/37) 2 ^(-5-3f/ + /) 

for the minus (—) sign in the last term in (1271) . and for U ^ —2, 

(-5 -3U + /)(-3 - 3E/ - v^/) 2 ^ 
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fc,E 4 



(C7 + 2) 2 (1 + 3f7 + /) 



(29) 



for the plus(+) sign. The above results reduce to the charged black hole case in the limit 
a — > 0, and ki, k 2 are functions of the angle 9 to be determined from appropriate boundary 
conditions. The above two solutions correspond to the so-called regular and accelerated 
branches pUl [22]. Once U is obtained, integrating (1251) using Z = dP/dr we get the ex- 
pression for P(r, 9) given below. Then the scalar function H(r, 9) for the metric is given via 

G22D- 

P(r, 9) = - - (r 4 + 6a 2 r 2 cos 2 9 - 3a 4 cos 4 9) - |m 2 / dr r 2 (l + a) 2 U(r, 9) (30) 

T T L J 

For the two special solutions, U — and U = —2, they correspond to A = and A = m 2 in 
flHUl) and have Ricci scalar values R = and R = 12m 2 , respectively. This coincides with 
the values we obtained in the previous section for spacetime with constant scalar curvature. 



Next, in order to check the physical properties of the metric given by ([IB"]) we want 
to transform the Kerr-Schild form to the Boyer-Lindquist coordinates. For the case with 
constant curvature, the general Boyer-Lindquist type transformation is given by [13] 

dT = dt+ (l-\r*)A r dr 

cltP 

dip = dip — Xadt + — - dr, (31) 

(H + a z )A r 

where A r (r, 9) = (r 2 + a 2 )(l — Ar 2 ) — rP(r, 9). The Maxwell equation (114p remains satisfied 
in the transformed Boyer-Lindquist coordinates, if the potential one-form is given as follows. 



AfM dx» = ^ (-dt + ^^-d(f> ) (32) 



Under the transformation fl3Tl) . the metric (|18[) takes the form 

, 9 A, r /, asin 2 ^ ,,\ 2 S , 9 S A e sin 2 9 , , 

ds 2 = --f — d0 + — dr 2 + —d9 2 + — °— adt - ± — - — J -d<p . (33) 



adt — - — ^ — -rfi 



Now the constraint equation (|2"3"1) for the function P(r, 9) is not preserved with the trans- 
formed metric ( 133]) if P has the #-angle dependence. Since we want to keep the solutions 
obtained in the Kerr-Schild form valid even in the transformed Boyer-Lindquist coordinates, 
we want to restrict ourselves to the cases of ^-independent P. From (1HD1 . we see that P 
becomes independent of 9 only in the following two cases, A = 0, U = and A = m 2 , U = —2, 
and the both give the same scalar function H, 

HM-VH^L, (34) 
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where M is the total mass of the black hole, and boundary conditions are assumed to be 
implemented. 

Here, we note that the metric (I33p with the scalar function (I34p and A = m 2 yields an 
exact Kerr-Newman-de Sitter solution with the cosmo logical constant A = 3m;?, where the 
crossover scale r c = m~ x behaves as an effective cosmological radius. 

The self-accelerated branch of the Schwarzschild solution in the DGP model [TUJ [22] 
exhibits the negative mass behavior, as its asymptotic behavior P/r « —f 2 M /r 2 + m 2 c r 2 at 
the scale r ^> r* shows. The term —f 2 M /r 2 looks like a 5D negative mass in the 5D viewpoint, 
which suggests that there may be some non-perturbative instabilities in this branch. In our 
constant curvature de-Sitter solution, the metric solution does not show such negative mass 
behavior at any scale as we see below. The Newton potential $ which corresponds to —P/2r 
in the Schwarzschild case above can be obtained from the metric solution fl33l) as 



-2$ = - - — (2Mr-Q 2 )+ m 2 c (r 2 + a 2 sin 2 9). (35) 

r z + cr cos^ a 

The first term is the usual positive mass contribution to the potential, the second term is the 
usual charge contribution, and no negative mass contribution appears in the above potential. 
Finally, in the de-Sitter solution the governing equation for horizon radius is given by 

A r = (r 2 + a 2 )(l - m 2 r 2 ) - 2Mr + Q 2 = 0. (36) 

and roHi provided that the total mass M lies in the range M le < M < M 2e where M le and 
M 2e are given by 



Mie = r yl + m 2 (33a 2 + 36Q 2 )(1 - m 2 c a 2 ) - (1 - Um 2 c a 2 - 12m 2 Q 2 + m^a 4 ) 3 / 2 , 

M 2e = —= — sj\ + m 2 (33a 2 + 36g 2 )(l - m 2 a 2 ) + (1 - 14m 2 a 2 - 12m 2 Q 2 + m 4 a 4 ) 3 / 2 . 

(37) 

Here tch, which is smaller than the crossover scale r c , is a cosmological horizon, and r + , r_ 
are outer and inner horizons, respectively. Note that the horizons r± and tch always have 
real positive values if the total mass lies between the masses M le and M 2e . 



4 Cosmological solution with de-Sitter background 

In this section, we will explicitly show that in the non-rotating limit our de-Sitter solution 
matches with a cosmological solution with accelerated expansion. In showing this, we follow 
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the method of matching conditions used in Ref. [15]. We shall see that the solution with 
A = m 2 c as the accelerating universe with the Hubble parameter 7i = r^ 1 in the non-rotating 
limit. In [T3], requiring that the underlying gravity theory respect Birkhoff's law the modified 
gravitational force law needed to generate the given cosmology was derived. Here we reverse 
the process: We examine the cosmological evolution corresponding to our metric (l34"l) with 
A = m 2 assuming the Birkhoff's law in the non-rotating limit. 

Following the approach of Ref. [T3] we consider a uniform sphere of dust and radia- 
tion(charge) and imagine that the evolution inside the sphere is exactly cosmological, while 
outside the sphere is filled with nonzero electric field. The mass and charge as the source are 
assumed to be unchanged throughout its time evolution. We take the cosmological metric 
to be the Robertson- Walker metric (flat space) inside the sphere: 

ds 2 = —dt 2 + a 2 {t)5ijdx l dx^ 

= -dt 2 + a 2 (t) [dr] 2 + rfdtt 2 } , (38) 

where dVL 2 is two-sphere. Outside the sphere we take a Schwarzschild-like metric 

ds 2 = -g 00 (r)dT 2 + g„{r)dr 2 + r 2 dtt 2 , (39) 

and rewrite it in a new form 

ds 2 = -N 2 (t, r])dt 2 + a 2 (t)dr] 2 + rW, (40) 

where r = r(t,T]). Now we set the coordinate transformation between (1391) and (1401) by 

dT = Tdt + T'drj, dr = rdt + r'dr], (41) 

where the dot and prime denote partial differentiations with respect to t and T], respectively. 
Rewrite ( 1391) under the transformation (T4T]) and compare it with ( T40i) . we obtain the following 
set of relations: 

N 2 = g 00 f 2 - g rr r 2 , 

= g 00 TT' - g„rr' ', 
-a 2 = gooT' 2 - g rr r' 2 . (42) 

In order to match the metric (l38l) inside the sphere with the metric (f4*0l) outside, we take the 
bounding surface (boundary) of the sphere to be fixed at rj = 1]*. By requiring continuity 
and smoothness at the boundary, we obtain the following matching boundary conditions: 

r(t,r]*) = i]*a(t), 
r'(t,rj # ) = a(t), 
N(t, V *) = 1, 

N'(t,ri*)=0. (43) 
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Applying the above matching conditions to (H2|) . the components of the metric in (1391) are 
determined in terms of a(t): 

g 00 = C 2 {1- tfa 2 ), g- 1 = 1 - r,l<x\ (44) 

where C 2 is an integration constant. 

In the non-rotating limit, from (1331) with ( 1341) for A = m 2 the components of our metric 
outside the sphere can be written as 

-1 2M 2 2 Q 2 

900 = 9rr = 1 m c r + ( 45 ) 

Comparing (jHJ) with ( T4"5l) and setting C 2 = 1, we obtain a Friedmann-like equation for the 
cosmological scale factor a(t): 

H 2 = (^) 2 = 2MCV 3 + m\ - Q 2 V ;'a-\ (46) 

When the scale factor a(t) goes to infinity at large time, the Hubble parameter 7i approaches 
the limit value Tio = m c = r~ x and yields the de Sitter expansion: 

a(t) oc e Hot . (47) 

Note that the cosmological metric (1381) with (1471) inside the sphere yields a constant Ricci 
scalar 12m 2 at large time, which is the same value from our metric 



5 Conclusion 

In this paper we investigate the Kerr-Newman solution on the braneworld with de-Sitter 
background in the DGP model . Assuming a ^-symmetry across the brane and a sta- 
tionary and axisymmetric metric ansatz in the Kerr-Schild form with a constant curvature 
background on the brane, we solve the constraint equation from (4+l)-dimensional gravity. 

When we restrict the spacetime with constant curvature, only flat and a de-Sitter cases 
are allowed by the constraint equation. When we transform the metric obtained in the 
Kerr-Schild form to the Boyer-Lindquist coordinates, the constraint equation in terms of the 
metric scalar function is preserved (so is the solution) only when the scalar curvature has 
two specific values, and 12m 2 . This is in accord with the restriction from the constraint 
equation since a de-Sitter type metric is used as the background metric in the Kerr-Schild 
form. 
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Unlike the self-accelerated branch of the Schwarzschild solution in the DGP model, our 
de-Sitter type solution does not have a negative mass behavior. This is checked by the 
inspection of our metric solution. About the effect of the bulk on the brane solution, we 
checked that there is no gravitational effect from the extra dimension in our constant curva- 
ture de-Sitter solution. We show this by evaluating the projected Weyl tensor(on the brane) 
which turns out to be only proportional to the energy-momentum tensor on the brane. 

In conclusion, we find an exact solution of Kerr-Newman type in a de-Sitter background 
in the Kerr-Schild form, which is also well behaved under the Boyer-Lindquist transforma- 
tion. Using the method of matching boundary conditions in the non-rotating limit, we show 
that this solution has the characteristic of cosmological solution with accelerated expansion. 



Acknowledgments 

The authors thank KIAS for hospitality during the time that this work was done. D.L. 
was supported by the Astrophysical Research Center of the Structure and Evolution of the 
Cosmos (ARCSEC), which is supported by the Korea Science and Engineering Foundation 
(KOSEF). 

References 

[1] S. Perlmutter et al, Astrophys. J. 517 (1999) 565. 

[2] L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 3370; 83 (1999) 4690. 

[3] G. Dvali, G. Gabadadze, and M. Porrati, Phys. Lett. B485 (2000) 208. 

[4] C. Deffayet, Phys. Lett. B502 (2001) 199; C. Deffayet, G. Dvali, and G. Gabadadze, 
Phys. Rev. D65 (2002) 044023. 

[5] A. Gruginov, New. Astron. Rev. 10 (2005) 311. 

[6] M. Porrati, Phys. Lett. B534 (2002) 209. 

[7] A. Lue and G. Starkman, Phys. Rev. D67 (2003) 064002. 

[8] C. Middleton and G. Siopsis, Mod. Phys. Lett. A19 (2004) 2259. 



12 



[9] A. Nicolis and R. Rattazzi, JHEP 0406 (2004) 059. 

[10] G. Gabadadze and A. Iglesias, Phys. Rev. D72 (2005) 084024. 

[11] Ee C.-Y. and D. Lee, Phys. Lett. B659 (2008) 58. 

[12] D. Lee, Ee C.-Y., and M. Yoon. la7Xiv":0711.1074 [hep-th]. 

[13] R. P. Kerr and A. Schild, Proc. Symp. Appl. Math 17 (1965) 199. 

[14] G. W. Gibbons, H. Lu, D. N. Page, and C. N. Pope, J. Geom. Phys. 53, 49 (2005). 

[15] A. Lue, R. S. Scoccimarro, and G. Starkman, Phys. Rev. D69 (2004) 044005. 

[16] A. N. Aliev and A. E. Gumrukcuoglu, Class. Quantum Grav. 21 (2004) 5081. 

[17] C. Deffayet, G. Gabadadze, and A. Iglesias, JCAP 0608 (2006) 012. 

[18] W. Israel, Nuovo. Cim. B44 (1966) 1; B48 (1967) 463(E). 

[19] T. Shiromizu, K. Maeda, and M. Sasaki, Phys. Rev. D62 (2000) 024012. 

[20] A. N. Aliev and A. E. Gumrukcuoglu, Phys. Rev. D71 (2005) 104027. 

[21] T. Dereli and M. Gurses, Phys. Lett. B171 (1986) 209. 

[22] G. Gabadadze and A. Iglesias, Phys. Lett. B632 (2006) 617. 

[23] M. H. Dehghani, Phys. Rev. D65 (2002) 104003. 



13 



